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ABSTRACT  (UNCLASSIFIED) 


A  novel  method  of  array  antenna  pattern  synthesis  is  presented  that  is  based  on  an  optimization  technique. 
The  complex  excitation  values  of  the  current  distribution  of  the  array  elements  are  taken  as  parameters.  The 
radiated  far  field  pattern  is  matched  to  one  that  is  prescribed  by  minimizing  a  global  mean  square  error.  The 
method  is  suitable  for  general  arrays  like  linear,  planar  and  conformal  ones,  and  does  not  require  the 
factorization  into  an  element  factor  and  an  array  factor.  In  contrast  with  Fourier-like  methods  the  array 
elements  may  be  arbitrarily  placed.  Although  the  theory  is  developed  for  arrays  of  identical  elements,  it 
applies  to  the  case  of  arrays  containing  non-identical  elements  as  well.  Numerical  results  are  presented  for 
unequally-spaced  planar  arrays  of  short  electric  dipoles,  x 
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SAMENVATTING  (ONGERUBRICEERD) 

Een  nieuwe  methode  voor  array  antenne  patroonsynthese  wordt  gepresenteerd,  die  gebaseerd  is  op  een 
optimaliseringstechniek.  De  complexe  excitatie-coefficienten  van  de  array  elementen  zijn  de  ontwerp 
parameters.  Het  gewenste  verre  veld  antennepatroon  wordt  zo  goed  mogelijk  benaderd  op  basis  het 
minimaliseren  van  een  gemiddelde  kwadratische  font. 

De  methode  is  ge.schikt  voor  zowel  lineaire,  vlakke  als  niet-vlakke  arrays  en  kent  geen  uitsplitsing  van  het 
verre  veld  patroon  in  een  element  en  array  factor.  In  tegenstelling  tot  Fourier-achtige  methoden  mogen  de 
antenne-elementen  willekeurig  geplaatst  worden.  Hoewel  de  theorie  is  opgezet  voor  arrays  van  identieke 
elementen  is  de  theorie  ook  bruikbaar  voor  niet-identieke  elementen.  Numerieke  resultaten  worden  getoond 
voor  een  vlakke  array  antenne  van  niet-equidistante  elektrische  dipolen. 
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1  INTRODUCTION 

In  array  antenna  synthesis  problems  an  optimum  choice  of  the  design  parameters  is  aimed  at  such 
that  the  radiation  pattern  of  the  array  approximates  in  a  certain  sense  a  desired  one.  As  the  design 
parameters  for  an  array  antenna  there  are  available  the  number  of  elements,  their  orientations  and 
positions  and  the  amplitudes  and  phases  of  their  excitations.  The  desired  radiation  pattern  is  either 
partly  fixed  by  one  or  more  of  its  characteristics  like  the  directive  gain,  the  beamwidth,  the 
sidelobe  level  or  directions  in  which  zero  power  has  to  be  radiated  (null  directions)  [l]-[6],  or 
specified  in  all  observational  directions  by  prescribing  a  complete  radiation  pattern.  In  the  latter 
case  we  are  dealing  with  pattern  matching. 

The  different  methods  of  pattern  matching  can  be  distinguished  by  the  manner  in  which  they 
model  the  synthesis  problem.  First,  there  are  the  sampling  methods  where  the  design  parameters 
of  the  array  antenna  are  chosen  to  match  (exactly  or  within  some  tolerance)  the  prescribed  pattern 
in  a  discrete  number  of  directions  of  radiation  [7]-[9].  Next,  there  are  many  methods  that  model 
the  matching  problem  as  a  minimization  problem.  In  this  case  the  design  parameters  are  chosen  to 
minimize  some  error  within  the  concept  of  a  given  error  criterion.  Commonly  used  error  criteria 
in  this  respect  are  the  minimax  error  [10]-[13]  and  the  mean  square  error  [4],[14]-[16].  Recently, 
methods  have  been  developed  that  model  the  matching  problem  as  an  intersection  problem  in  the 
space  of  admissible  radiation  patterns  [17]-[19].  Each  requirement  of  the  radiation  pattern  that  is 
to  be  synthesized  leads  to  the  admissibility  of  a  subset  in  the  set  of  all  the  realizable  radiation 
patterns.  A  radiation  pattern  that  is  an  element  of  the  intersection  of  all  these  subsets  meetr  all  the 
requirements  stated  and  is  looked  for. 

Many  of  the  array  synthesis  methods  apply  to  equispaced  linear  or  planar  arays.  However, 
unequally-spaced  arrays  are  also  of  great  interest  [13]-[15],  [20]-[27].  In  case  the  array  is  to  be 
used  as  a  directional  antenna  with  high  resolution  rather  than  high  gain,  unequally-spaced  arrays 
will  often  be  superior  to  equispaced  arrays.  The  unequally-spaced  arra;^  s  require  fewer  elements 
to  produce  a  certain  resolution,  and  their  main  beam  can  be  steLred  over  a  larger  frequency 
bandwidth  than  is  the  case  with  an  equispaced  array.  Arrays  w’tn  variable  interelement  spacings 
can  approximate  a  radiation  pattern  more  closely  than  arrays  with  constant  spacings  because  of 
the  extra  degree  of  freedom  that  is  introduced  by  the  variability  of  the  spacings. 

In  most  methods  known  for  unequally-spaced  arrays  the  complex  excitations  of  the  elements  of 
the  array  are  fixed  and  their  relative  positions  ar^  taken  as  the  design  parameters  [20]-[27].  Once 
these  positions  have  been  fixed  it  may  be  advantageous  to  perform  a  second  step  in  the  synthesis 
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procedure  by  considering  the  complex  excitations  as  the  design  parameters  and  apply  a  method  to 
find  the  optimum  values  of  the  latter  in  order  to  further  improve  the  ^proximation  to  the 
prescribed  radiation  pattern.  So,  in  this  case  we  need  a  synthesis  method  that  performs  a  pattern 
matching  given  the  positions  of  the  fixed,  unequally-spaced  array  elements.  For  unequally-spaced 
linear  arrays  such  a  method  has  been  described  in  [13].  For  unequally-spaced  planar  arrays  the 
iterative  sampling  method  [8,9]  or  the  computational  methods  described  in  [15]  are  suitable.  A 
disadvantage  of  the  iterative  sampling  method  is  that  there  is  no  guarantee  of  the  existence  of  a 
solution.  In  case  a  solution  is  found,  this  need  not  be  an  optimum  one  in  any  sense.  Thus,  better 
solutions  may  exist  that  are  not  found  by  this  method.  The  method  described  in  [15]  does  not  have 
these  difficulties  since  it  is  based  on  the  minimization  of  an  error  function.  It  provides  an  unique 
optimum  solution.  However,  that  theory  has  only  been  applied  to  point  soinces  that  are  arbitrarily 
distributed  in  a  plane  and  to  a  pattern  synthesis  within  that  plane.  In  addition,  the  error  function  is 
based  on  the  deviations  of  the  synthesized  radiation  pattern  from  the  prescribed  one  at  a  finite  set 
of  discrete  points  of  observation.  There  is  no  control  of  the  radiation  field  at  the  points  in 
between. 

In  this  paper,  a  pattern  synthesis  method  is  presented  for  unequally-spaced  array  antennas  that 
may  consist  of  non-identical  elements.  It  involves  the  minimization  of  the  quadratic  error 
integrated  over  the  unit  sphere  on  which  the  far-field  radiation  pattern  is  defined,  and  applies  to 
the  (transverse)  electric  field  radiation  characteristic  defined  on  that  sphere.  Since  the  (transverse) 
magnetic  field  radiation  characteristic  is  identical  to  the  one  for  the  electric  field,  the  mismatch  in 
the  magnetic  field  is  simultaneously  minimized  with  the  mismatch  in  the  electric  field.  The  class 
of  admissible  patterns  that  may  be  prescribed  is  chosen  in  accordance  with  the  physics  of  the 
electromagnetic  radiation  of  the  array  elements. 

The  analysis  is  carried  out  in  the  frequency  domain  (co-domain)  and  considers  arrays  consisting  of 
arbitrarily  oriented,  possibly  different,  elements  (electric  dipoles,  radiating  apertures).  The  action 
of  each  element  is  characterized  by  a  properly  normalized  spatial  distribution  function  of  volume 
or  surface  current,  and  a  complex  excitation  coefficient  that  determines  the  amplitude  and  phase 
impressed  on  this  distribution.  The  excitation  coefficients  are  the  "design  parameters"  of  the 
synthesis  problem.  It  is  assumed  that  they  can  be  chosen  independently,  i.e.,  mutual  coupling  is 
neglected.  The  radiation  field  of  the  array  antenna  follows  by  the  summation  of  the  fields 
generated  by  the  individual  elements.  In  this  representation,  the  complex  excitation  coefficients 
occur  as  expansion  coefficients.  Next,  the  representation  is  substituted  into  the  error  criterion  and 
a  variational  method  is  applied  to  minimize  it.  This  leads  to  a  system  of  equations  the  solution  of 
which  jjTovides  the  optimum  choices  of  the  complex  excitation  coefficients. 
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The  theory  is  more  general  than  for  the  case  of  identical  elements  which  is  usually  considered  and 
where  a  factorization  into  an  element  factor  and  an  array  factor  is  carried  through.  The  more 
general  case  that  we  consider  admits  no  such  factorization.  (Note  that  the  removal  of  this 
restriction  yields  an  extra  degree  of  freedom  in  the  design  of  the  antenna.) 

The  array  antenna  is  assumed  to  be  present  in  free  space.  Numerical  results  are  presented  for 
different  planar  arrays.  The  analysis  is  organized  such  that  a  number  of  steps  can  be  carried  out 
analytically,  which  leads  to  a  saving  in  computation  dme. 
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2  BASIC  EQUATIONS  AND  SOURCE-TYPE  REPRESENTATIONS  FOR  THE 

ELECTRIC  AND  THE  MAGNETIC  FIELD  STRENGTHS 

The  array  antenna  is  present  in  free  space  having  electromagnetic  medium  properties  given  by  the 
permittivity  Eq  and  the  permeability  To  specify  the  position  in  the  configuration  we  employ 
the  standard  coordinates  respect  to  a  fixed,  orthogonal,  Cartesian  reference  frame 

with  origin  O  (the  phase  center  of  the  array  antenna)  and  the  three  mutually  perpendicular  base 
vectors  of  unit  length  each.  In  the  indicated  order  the  base  vectors  form  a  right-handed 

system.  The  subscript  notation  for  Cartesian  vectors  and  tensors  is  employed  and  the  summation 
convention  applies.  The  corresponding  lower-case  subscripts  are  to  be  assigned  the  values  1,2,3. 
This  usage  leads  to  expressions  that  are  effortlessly  translated  into  the  statements  of  a  computer 
program  written  in  some  high-level  programming  language.  Partial  differentiation  is  denoted  by 
3;  dp  denotes  differentiation  with  respect  to  Xp. 

The  analysis  is  performed  in  the  frequency  domain  (o>-domain).  The  complex  time  factor  is 
expOoJt).  where  t  is  the  time  coordinate  and  p^-\.  The  electromagnetic  wave  field  generated  by 
the  array  antenna  is  represented  by  the  electric  field  strength  and  the  magnetic  field  strength  //, 
and  the  volume  source  densities  by  and  Kj.  The  corresponding  Maxwell  equations  are  [28,  p.3] 


(1) 

(2) 

where  Ej^^p  is  the  completely  anti-symmetric  unit  tensor  of  rank  three  (Levi-Civita  tensor)  defined 
by 


^kmp 


-1  when  {k.m,p]  is  an  odd  permutation  of  {1,2,3}, 
0  when  not  all  subscripts  are  different, 

1  when  { k,m,p  1  is  an  even  permutation  of  { 1 ,2,3  j . 


(3) 


Equations  (1)  and  (2)  govern  the  propagation  of  the  electromagnetic  waves  in  free  space 
generated  by  the  source  distributions  7^  and  Ky 
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The  frequency-domain  source-type  representations  for  the  electric  and  the  magnetic  field 
strengths  of  the  field  radiated  by  distributed  sources  present  in  the  bounded  source  domain  are 
given  by  ([28,  p.38]) 

Ek  =-;COMi+  W 

and 

=-j(SX^j+  +  (5) 

in  which 

r 

Ai(Ar,co)=  I  ,  G(jc-j:  ,(0>/*(jc  ,CD)dV  (6) 

JjceD^ 

is  the  electric  source  current  vector  potential  and 

G(jr-x',(o)/s:,{Jc’,co)dV  (7) 


Fj(x,w)  = 


is  the  magnetic  source  current  vector  potential.  Here, 

/7(x,co)  =  expC-yZ'oIxIV^Jtlxl,  (8) 

is  the  scalar  free-space  Green's  function  (point-source  solution)  with 

ko  =  co(po£o)  =  to/Cfl  =  2nAo,  (9) 

where  Cq  is  the  electromagnetic  wavespeed  in  vacuum  and  Xq  is  the  free  space  wavelength.  In  our 
pattern  synthesis  problem  we  are  only  interested  in  the  radiation  field  in  the  far-field  region  for 
which  the  far-field  parts  provide  the  relevant  expressions.  In  these  expressions  the  following 
approximations  are  used;  (Collin  and  Zucker  1969,  p.l40) 

ix-xi  =  Ixl-^^/  +  vanishing  terms  as  1x1— >«  (9a) 


in  the  exponential  function  occurring  in  the  Green’s  function,  where 
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=  jfj/lxl  (9b) 

is  the  unit  vector  in  the  direction  of  observation,  and 

Ix-x'l  =  Ixl,  (9c) 

in  each  denominator  of  the  expression  for  the  far-field  part.  The  far-field  approximation  follows 
as  ([28],  p.39]) 

{ )  (JT.CO)  =  { ]  (|,(o)exp(-;*o'^')/47ibrl  as  ( 10) 

in  which  ^=x/lxl  is  the  unit  vector  in  the  direction  of  observation, 

£r=  (H) 

and 

W  “  =  -jmo{^jp-^ji,p)Fp-  -  ( 1 2) 

where 


(/l,‘”/'p~](|,a))  = 


xeD^ 


(13) 


are  the  far-field  radiation  characteristics  of  the  current  distribution. 

As  (10)  shows,  the  electric  and  the  magnetic  field  strength  have,  in  the  far-field  region,  the 
structure  of  a  spherical  wave  that  expands  radially  from  the  origin  of  the  coordinate  system 
(which  is  also  denoted  as  the  phase  center  of  the  far-field  approximation),  the  latter  being  chosen 
in  the  neighbourhood  of  the  radiating  sources,  with  an  amplimde  that  depends  on  the  direction  of 
observation  and  that  decreases  inversely  proportional  to  the  distance  from  the  sources.  The 
amplitude  radiation  characteristics  depend  only  on  the  direction  of  observation  and 

on  to. 

The  far-field  amplitude  radiation  characteristics  E/^"  and  //y-  of  the  electric  and  the  magnetic  field 
strengths  are  not  independent  of  each  other.  It  is  easily  verified  that  the  right-hand  sides  of  (11) 
and  (12)  are  interrelated  in  the  following  way:  (De  Hoop  1990,  p  10-27) 
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^kmpznfip°°+  T£r  =  0, 

(14) 

-e;m4m£“+  ZHf  =  0, 

(15) 

in  which 


Y  =  (€o/iio)*^  and  Z  =  (iVeo)*'^ 


(16) 


are  the  plane-wave  admittance  and  impedance  of  free  space,  respectively.  Note  that 


4^r  =  0and^//^“  =  0,  (17) 

i.e.,  in  the  far-field  region  the  electric  and  the  magnetic  field  strengths  are  transverse  with  respect 
to  the  local  radial  direction  of  propagation  ^  of  the  wave,  and  are,  since  6^  and  have  unit 
magnitudes,  proportional  with  proportionality  factors  Y  and  Z.  The  directive  gain  (or  directivity) 
of  tlie  antenna  array  [28,  p.33]  is 


D(|.co) 


47t£'y~(|,(0)£/-(|,(0) 

~‘A  ■  - - 

£,~(|',(0)£/-(|',co)dQ 


4n//,~(|,to)///-(|,(D) 
7 - 


for 


(18) 


where  is  the  unit  sphere  around  the  phase  center  of  the  antenna. 
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3  OPTIMIZATION  TECHNIQUES  IN  ARRAY  ANTENNA  PATTERN 

SYNTHESIS 

3. 1  Introduction 

In  this  section,  the  optimization  technique  will  be  developed  for  elements  with  electr'  current 
volume  source  distributions  and  for  a  prescribed  radiation  characteristic  of  the  electric  field 
strength  within  the  class  of  vector  functions  that  are  tangential  to  the  unit  sphere.  The  technique 
for  other  types  of  elements  (such  as  radiating  apertures)  goes  along  the  same  hnes,  while  in  view 
of  the  equahty  of  the  (transverse)  electric  and  magnetic  field  radiation  characteristics,  the 
minimization  of  the  mismatch  in  electric  field  simultaneously  minimizes  the  mismatch  of  the 
corresponding  magnetic  field.  The  volume  source  current  differs  from  zero  in  a  given  bounded 
domain  is  the  union  of  the  domains  {f>^(/);  /=!,.. .A^}  occupied  by  the  separate  elements 

(/V  =  number  of  elements). 

3.2  The  Pattern  Synthesis  Problem  formulated  as  a  Minimization  Problem 

To  arrive  at  the  error  criterion  that  will  be  minimized  in  our  synthesis  problem,  we  first  introduce 
the  inner  product  <.,.>n  defined  on  the  unit  sphere  Q.  For  two  vector  functions  R(^,co)  and  0(^,co) 
this  inner  product  is  taken  as 


<PQ>n  = 


a(|)P^(|,co)ep-(|.®)dQ, 


(19) 


where  is  the  unit  sphere  around  the  origin  and  w(^)  is  a  real  positive  weighting  function 
(h  (^)>0  for  all  ^)  and  *  denotes  complex  conjugate.  The  global  error  function  that  describes  the 
mismatch  between  the  prescribed  radiation  characteristic  and  the  radiation 

characteristic  £*  realized  by  the  array  is  chosen  as 


ERR(£")  =  ![n  =  (<£“’-£D'-J?'><>-£0~>n)//2. 


(20) 
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The  weighting  function  occurring  in  the  definition  of  this  norm  (cf.  19)  provides  the  possibility  to 
match  the  radiation  pattern  of  the  antenna  in  some  directional  area  better  than  in  other  directional 
areas. 

For  the  class  of  elements  under  consideration,  we  have  in  view  of  (1 1) 


where  is  given  by  (cf.  (13)): 


(21) 


A,“(|,co)  = 


r 


xeD^ 


J^(x,(o)exp(Jk(^^,)dV. 


(22) 


Note  that  the  integration  variable  x  has  been  replaced  by  x. 

Let  j(x,wj)  be  the  normalized  source  distribution  of  the  element  occupying  the  domain 
(/=!,...//),  and  let  c(O)',/)  be  its  complex  excitation  coefficient,  then 


N 


(23) 


with 


a,.“(^.to;/)  = 


xeom) 

for  7=1,. ..,N. 


;V(j:,co;/)exp(/'7:o^,ar^)dF, 


(24) 


and 


j  c(co;/)et~(|.(otO 


(25) 


with 
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=  -/•co^o(5tr-4/fct)ar“(4.“;^  (26) 

for/=l,...//. 

Substitution  of  the  expansion  (25)  into  (20)  gives  for  the  squared  norm  of  the  mismatch  in  the 
field 


\E--ED- 11 2  „=  |j£D~  Ij2  ^ .  2Re[£  ^c*(co;/)<E^“(|,to),e“(|,co;/)>] 
N  N 

+  I  ^  I _  c(co;/)c*(to;/)<e~(|,tD;/),e“(|,co;y)>. 


(27) 


To  arrive  at  the  condition  for  the  optimum  complex  excitation  coefficients  c®P'(o>;/),  that 

minimize  the  squared  error  (27),  we  apply  a  variational  analysis.  As  a  result,  the  following  system 
of  equations  is  obtained  (see  Appendix  A): 


N 

L  cop‘(co;/)<c“’(^,co;/),e“(§,<o;y)>  =  <£:®“(^,C0),e~(^,(0;7)> 

for/=l,...A^.  (28) 

Given  an  array  antenna  for  which  the  electric  current  volume  source  distributions  of  its  elements, 
and  hence  the  functions  jix,(Oj),  are  known,  the  matrix  elements  of  the  system  matrix  in 

(28)  follow  from  (24)  and  (26).  Together  with  the  specification  of  the  right-hand  side  through  the 
desired  pattern  the  system  of  equations  (28)  is  completely  known  and  can  be  solved  for  the 
optimum  excitations  c^P'(co;/),  The  prescribed  radiation  characteristic  of  the  electric  field 

strength  is  taken  from  the  class  of  vector  functions  that  are  tangential  to  the  unit  sphere  (cf.  (21)), 
i.e.. 


(29) 

where  A,.®”  is  user  prescribed.  In  fact,  we  prescribe  the  far-field  radiation  characteristic  of  the 
electric  current  distribution  A^^”  from  which  follows  that  is  to  be  approximated  by  E^°°. 
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3.3  Construction  of  The  System  of  Equations 

The  elements  of  the  system  matrix  m  (28)  are,  with  the  aid  of  (19),  (24)  and  (26),  found  as 


<e“(|,w;/),  e“(|,co;7)> 


xeDT{r) 


jp(x,OiJ) 


fOT 


where 


(30) 


Fp^(xj,(si)  = 


£6  Q  - 


(31) 


and 


Wp^(^,(M)  =  (COPo)2w(|)(5p<,-£^).  (32) 

The  components  of  the  right-hand  side  vector  in  (28)  have,  for  the  class  of  functions  (29),  the 
form 


<£^~(£,(0),e~(£,co;y)> 


for 


(33) 


where 


C,C>’.a))  = 


^p<?(|'®)Ap^~(|-®)exp(->*:o^^^)d£2. 


(34) 


Through  (30)  -  (34)  and  (28)  the  optimization  problem  can  be  solved.  In  the  following  section 
several  test  cases,  when  the  elements  are  electric  dipoles,  will  be  considered. 
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4  TEST  ARRAYS 

Our  pattern  synthesis  method  is  illustrated  for  a  number  of  planar  arrays  of  ND  short  electric 
dipoles.  The  dipoles  are  taken  to  be  uniformly  oriented  and  non-uniformly  placed  in  the  Xj^2~ 
plane.  The  position  vector  of  the  barycenter  of  the  /-th  dipole  element  is  denoted  by  xQ), 
1=1,.. and  the  uniform  orientation  of  the  electric  dipoles  by  the  unit  vector  a.  The  electric 
current  source  distribution  of  the  /-th  element  is  taken  as  [28,  p.29] 


j{x,(a,r)  =  fl6[jc-jc(/)]  for  I=l,...J^D,  (35) 

where  6(x)  is  the  three-dimensional  Dirac  delta  function.  Substitution  of  the  somce  distributions 
(35)  into  (30)  and  (33)  leads  to  the  system  of  equations 


ND 

l^c°v\(Si-J)apa/p^{x{D,y{J),(ii\  =  a^G^\y{T),(si\  for  (36) 

We  now  take  the  weighting  function  w(^)=l.  For  this  case,  an  analytical  expression  of  the 
function  Fp^  in  (3 1)  can  be  found  (see  Appendix  B);  it  is  given  by 


Fp^(jrO',co)  =  47t(a)Po)2 


^pq 


hpq-3 


iXp-yp)  (x,-y,) 

~R  R 

(Xp-yp)  ix^-y^) 


R 


R 


sm{koR)/{koR) 


(kgRyHcosikoR)  -  smik(^)/(koR)], 


(37) 


in  which  /?=ljr-yl.  The  components  a,G,[y(y),0)],  J=\,...JND,  of  the  right-hand  side  of  the  system 
of  equations  (36)  becomes,  with  the  choice  H'(^)=1, 


a<,G,l>(/).(0]=  (COPo)2j^^^ 

exp[-7l:^^j^-,(/)]dn,  for  /=  1  ,...ND. 


(38) 


TNO  report 


Page 

17 


In  this  expression  is  the  prescribed  far-field  electric  source  current  vector  potential.  Because 
the  anteima  configuration  is  symmetric  with  respect  to  the  Xj,X2-^\siOS,  a  vector  potential  is 
prescribed  that  also  reflects  this  symmetry  and  we  specify  its  directional  dependence  in  the 
halfspace  x^>0  only.  The  prescribed  far-field  electric  source  current  vector  potential  of  the  whole 
antenna  system  is  taken  to  be  linearly  polarized  in  a  constant  real  direction  L  and  non-zero 
in  a  small  area  about  the  principal  direction.  It  is  taken  as  a  function  of  through  (see 
Fig.  2, Chapter  VI,  Numerical  Results) 


=  Lp4jW[^j-cos(5)]  for  CX^j,  (39) 

where  Hit)  is  the  Heaviside  unit  step  function  defined  as  Hit)  =  {0,1/2,!}  when  {/<0,f5=0,r>0}  and 
5<7t/2  is  a  positive  constant  that  determines  the  beamwidth  of  the  prescribed  radiation  pattern. 
Substitution  of  (39)  into  (38)  leads,  considering  the  symmetry  of  the  integrand  with  respect  to  the 
X;^2*pl2ne,  to 


a,C,0(/),(i)]  =  2((Oiio)2 


^5//[^5-cos(6)]exp[-;^^^,(/)]dI2, 


for/=l,...//D, 


(40) 


where  Q'  denotes  the  half  of  the  unit  sphere  Q  located  in  ^j>0.  Expression  (40)  can  be  determined 
analytically.  To  this  end,  we  rewrite  (40)  as 


agG^\yiJ),w]  =  ia^^  +  ko-^a^pd^p)Tm\yiJ),(si\> 

fory=l,...A^D,  (41) 

where  9^  (/.- 1 ,2)  denotes  partial  differentiation  with  respect  to  y*(7)  and 


INTl>(/).a)]  =  2(a)^(,)^ 


for/=l,...jVD. 


(42) 


The  analytical  evaluation  of  (42)  leads  to  (see  Appendix  C) 
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INT[y(/),co] 


471(03^0)^ 

- ^  QJjiQ),  for/=l,...^D. 

(A:ol>(y)l)2 


(43) 


in  which  /y  is  the  Bessel  function  of  the  first  kind  and  of  order  one  and  its  argument  is  given  by 
G=^ol3’(-^lsin(8).  Substitution  of  (43)  into  (41)  and  carrying  out  the  different  partial 
differentiations  gives  us  the  expression  for  the  right-hand  side 


afi^\y{f)M  =  47c(co|i(,)2sin25 


JiiQ) 


OpLp  -  sin26 


yq{J)yp(J) 

\y{f)\^ 


-  sin25 


JiiQ) 


dpLp  -  ^qLp 


yq(J)yp(J) 
\y(J)\^ 


02 

fory=i,...;;^D, 

in  which  J2  is  the  Bessel  function  of  the  first  kind  and  of  order  two. 


(44) 


With  the  expressions  (37)  and  (44)  the  system  of  equations  (36)  is  analytically  known.  To  fix  it, 
the  positions  x(/),  7=1, ..WD,  and  the  uniform  direction  a  of  the  electric  dipoles  have  to  be  given 
and  the  state  of  polarization  and  beamwidth  of  the  radiation  pattern  that  is  to  be  matched  have  to 
be  prescribed.  For  some  unequally-spaced  arrays  of  electric  dipole  antennas  numerical  results  will 
be  given  in  the  next  section. 
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5  NUMERICAL  RESULTS 

For  the  generation  and  solution  of  the  complex  system  of  equations  (36),  a  Fortran  77  code  has 
been  written  called  SYNPA  (=  SYNthesise  of  Planar  Array).  The  input  data  of  SYNPA  are  the 
positions  and  orientations  of  the  electric  dipoles  and  the  width  and  the  state  of  polarization  of  the 
prescribed  beam.  As  output  SYNPA  produces  several  print  and  plot  files  containing  the  chosen 
antenna  configuration  (plot),  the  computed  complex  excitations  (print),  the  global  error 

NERR=ERR/||£0“||  (45) 

(print),  where  ERR  is  defined  in  (20)  and  ||£^^~||  is  the  norm  of  the  described  radiation  pattern, 
and  the  corresponding  radiation  pattern  (plots).  From  the  synthesized  radiation  field  2D-plots  of 
the  H-,  D-,  and  £-plane  as  well  as  a  3D-plot  are  given.  Additional  Fortran  77  codes  have  been 
written  that  have  these  plot  files  as  input  and  produce  screen  and  hardcopy  output.  These  codes 
make  use  of  the  subroutines  from  the  UNIRAS  library  [29].  As  an  additional  option,  one  can 
carry  out  the  synthesis  for  digitized  complex  excitations.  This  option  is  of  interest  when  we  have 
to  deal  with  n-bits  amplifiers  and  m-bits  phase  shifters. 

For  the  evaluation  of  the  right-hand  side  of  the  system  of  equations  (36)  SYNPA  makes  use  of  the 
Bessel  functions  from  the  IMSL  Library  [30].  For  the  solution  of  (36)  SYNPA  uses  the  conjugate 
gradient  method  presented  in  [31]  without  any  preconditioning.  SYTJPA  has  been  implemented  on 
a  Vax  6310  computer  and  computation  times  needed  for  the  numerical  examples  given  in  this 
section  are  in  the  order  of  seconds. 

The  synthesis  problem  is  illustrated  for  a  number  of  square,  planar  arrays  of  9x9  electric  dipoles. 
For  each  array,  the  center  element  is  taken  to  be  at  the  origin  O  of  the  reference  frame  (phase 
center  of  the  antenna)  and  the  elements  are  arranged  such  that  they  are  parallel  to  the  coordinate 
axes.  The  currents  in  the  electric  dipoles  flow  parallel  to  the  Xy-direction.  The  (non-uniform) 
arrangements  of  the  elements  in  the  Xj-  and  X2-directions  are  chosen  to  be  identical  and  in  such  a 
way  that  the  array  configuration  is  symmetric  with  respect  to  the  two  coordinate  axes,  (see  Fig.  1). 
The  element  rows  parallel  to  the  xy-axis  are  numbered  as  -4,  -3,  ...,  3,  4,  where  row  0  is  the  one 
that  coincides  with  this  axis.  The  spacing  between  rows  i  and  j  is  denoted  by  and  is  expressed 
in  terms  of  wavelengths.  By  specifying  Sq/,  J;?,  S2s  and  Sj^,  the  element  spacings  of  the  whole 
array  are  fixed,  (see  Fig.  1).  The  numerical  experiments  are  carried  out  for  the  four  arrays  defined 
in  Table  I. 
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PLANAR  ARRAY  OF 
ELECTRIC  DIPOLE  ANTENNAS 

t 


-*■  -  Prescribed  far-field 


SCALE  IN  WAVELENGTHS 


Fig.  1 :  The  square  planar  array  antenna  of  9x9  unequally  spaced  electric  dipoles.  The  prescribed 

vector  potential  A  is  shown  from  which  the  prescribed  electric  field  E  follows. 


Results  for  grids  1-IV 

Grid 

^01 

Si  2 

S23 

I 

^4 

Type  of  grid 

Number  of 

NERR 

W 

(?t) 

W 

(A.) 

iterations 

(%) 

1 

0.5 

0.5 

0.5 

0.5 

equispaced 

23 

39 

II 

0.3 

0.3 

0.3 

0.3 

equispaced 

57 

43 

III 

0.5 

0.6 

0.7 

0.8 

unequally  spaced 

20 

37 

IV 

0.5 

0.75 

1.0 

1.25 

unequally  spaced 

19 

46 

Table  I:  The  four  different  grids  with  their  spacings  Sqj,  j  j2,  S23  ^34-  number  of  iterations 

needed  for  the  solution  of  the  corresponding  system  of  equations  for  each  grid  and  the 
related  computed  normalized  mean  square  error  NERR  arc  shown. 
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The  radiation  pattern  that  has  to  be  approximated  is  prescribed  via  the  far-field  electric  source 
current  distribution  Ar^°°  (cf.  (29),  see  Fig.  1).  We  choose  the  direction  of  this  vector  in  the 
positive  jTy-direclion.  Its  angular  dependence  is  given  by  (39)  in  which  we  choose  5  =  15”,  i.e.  the 
prescribed  beamwidth  is  30”.  In  Fig.  2  the  prescribed  radiation  field  is  presented. 


PRESCRIBED  |e| 

5M(e)cos(«) 


Prescribed  beam  width  =  30* 
Number  of  elements  =  8l 


Fig.  2;  The  prescribed  radiation  pattern  of  beamwidth  30”. 

The  results  of  the  synthesis  for  the  four  planar  arrays  are  given  in  Table  1  and  Figs.  3-10.  In  Table 
1  the  number  of  iteration  steps  needed  for  the  solution  of  the  system  of  equations  for  each  of  the 
four  grids  are  given.  A  large  number  of  iteration  steps  suggests  a  poor  condition  of  the  system 
matrix.  Further,  in  Table  I  the  normalized  global  mean  square  error  NERR  is  given  which  is  a 
measure  of  the  accuracy  of  approximation.  In  Figs.  3-6  three-dimensional  plots  of  the  synthesized 
radiation  fields  are  shown  and  in  Figs.  7-10  the  two-dimensional  plots  of  the  prescribed  and 
synthesized  radiation  fields  in  the  //-,  D-  and  £-planes  are  shown. 
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SYNTHESIZED  |Ei 

sKe)cos(«)  siN(e)siN(9i) 


Prescribed  beam  width  =  30* 
Number  of  elements  ^  81 


Fig.  3:  The  synthesized  total  radiation  pattern  corresponding  to  grid  I. 


SYNTHESIZED  |E| 

3f.4(8)COS(«)  SiN(©)siNC0) 


Prescribed  beam  width  *  30* 
Number  of  elements  «  81 


Fig.  4: 


The  synthesized  total  radiation  pattern  corresponding  to  grid  11. 
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SYNTHESIZED  |E| 

SM(0)cos(«)  siN<<9)sirNjC0) 


Prescribed  beam  width  -  30' 
Number  of  elements  «  61 


Fig.  5:  The  synthesized  total  radiation  pattern  corresponding  to  grid  III. 


SYNTHESIZED  IE| 

SM(e)cos(«)  siNCejsiNC^s) 


Prescribed  beam  width  e  30* 
Number  of  elements  e  81 


Fig.  6: 


The  synthesized  total  radiation  pattern  corresponding  to  grid  IV. 
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TOTAL  RADIATION  PATTERN 
IN  THE  H-.  D-  AND  E-PLANES 


—  =  PRESCRIBED  |  E  | 

H-,  D-  AND  E-PLANES 

SYNTHESIZED  |  E  |  IN 
H-PLANE 
--  D-PLANE 

—  E-PLANE 


Prescribed  beam  width  =  30“ 
Number  of  elements  =  81 


Fig.  7:  The  total  radiation  pattern  in  H-,  D-  and  £-planes  of  the  prescribed  and  synthesized  radiation 

field  corresponding  to  grid  I. 


TOTAL  RADIATION  PATTERN 
IN  THE  H-,  D-  AND  E-PLANES 


—  *  PRESCRIBED  |  E  | 

H-.  D-  AND  E-PLANES 

SYNTHESIZED  |  E  1  IN 
H-PLANE 

. D-PLANE 

- E-PLANE 


Number  of  elements  =  81 


Fig.  8: 


The  total  radiation  pattern  in  H-.  D-  and  £-planes  of  the  prescribed  and  synthesized  radiation 
field  corresponding  to  grid  II. 
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TOTAL  RADIATION  PATTERN 
IN  THE  H-.  D-  AND  E-PLANES 


—  =  PRESCRIBED  |  E  | 

H-.  D-  AND  E-PLANES 

SYNTHESIZED  |  E  |  IN 
H-PLANE 
. D-PLANE 

—  E-PLANE 


Prescribed  beam  width  =  30° 


Number  of  elements  =  81 


Fig.  9: 


The  total  radiation  pattern  in  H-,  D-  and  £-planes  of  the  prescribed  and  synthesized  radiation 
field  corresponding  to  grid  ID. 


LU 


TOTAL  RADIATION  PATTERN 
IN  THE  H-.  D-  AND  E-PLANES 


SIN(0) 


- PRESCRIBED  |  E  | 

H-,  D-  AND  E-PLANES 

SYNTHESIZED  |  E  |  IN 
H-PLANE 

-  -  -  D-PLANE 

—  E-PLANE 


Prescribed  beam  width  =  30° 
Number  of  elements  =  81 


Fig.  10; 


The  total  radiation  pattern  in  //-,  D-  and  f-planes  of  the  prescribed  and  synthesized  radiation 
field  corresponding  to  grid  IV. 
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Grids  I  and  II  are  in  fact  identical  grids  that  are  excited  at  different  frequencies.  For  grid  I  with  its 
0.5  wavelength  interelement  spacing  we  have  a  higher  frequency  of  excitation  than  for  grid  n 
with  its  0.3  wavelength  interelement  spacing.  From  Table  1  we  learn  that  the  global  error  NERR 
for  the  synthesis  for  grid  II  is  larger  than  the  one  for  grid  I.  Further,  the  plots  of  Figs.  3-4  and 
Figs.  7-8  show  that  the  sidelobes  of  the  synthesized  pattern  for  grid  n  are  higher  than  the  ones  for 
grid  I.  Finally,  the  number  of  iterations  needed  for  the  solution  of  the  system  of  equations  for  grid 
n  is  larger  than  the  one  for  grid  I.  Apparently  the  condition  of  the  system  matrix  in  ^be  former 
case  is  poorer. 

The  experiments  for  grid  HI  and  FV  are  the  ones  for  unequally-spaced  arrays.  The  results  for  grid 
ni  show  a  lower  global  error  NERR  than  the  ones  obtained  for  the  equally-spaced  grids  I  and  n. 
The  sidelobes  of  the  synthesized  radiation  pattern  for  grid  HI  are  higher  than  the  ones  tor  grid  I 
and  II  but  the  approximation  of  the  main  beam  is  better  than  for  the  grids  I  and  II  (see  Fig.  9).  In 
grid  rv  we  have  unequally-spaced  electric  dipoles  with  larger  interspacings  than  for  grid  EH.  For 
this  case  the  main  beam  is  even  better  approximated  (Fig.  6  and  Fig.  10).  However,  the  sidelobe 
level  is  much  higher  than  for  grid  IH  and  gives  rise  to  a  global  error  NERR  that  is  higher  than  for 
the  case  with  grid  HI. 

Although  Figs.  3-10  show  that  for  the  four  experiments  we  have  a  good  approximation  of  the 
prescribed  radiation  pattern  we  found  rather  large  global  errors  NERR  in  the  order  of  40%.  These 
large  errors  are  mainly  due  to  the  presence  of  the  steep  slopes  of  the  prescribed  beam.  A 
prescribed  beam  with  less  steep  slopes  and  the  same  beamwidth  will  result  into  a  significantly 
smaller  NERR.  However,  for  such  a  prescribed  radiation  pattern  the  right-hand  side  of  the  system 
of  equations  cannot  be  evaluated  analytically  and  a  numerical  evaluation  of  the  resulting  integrals 
is  required  (which  can  easily  be  carried  out). 

Next,  the  pattern  matching  is  carried  out  for  grid  I  (see  Table  1)  for  the  case  that  the  elements  are 
excited  with  2-bits  amplitude  excitations.  (The  number  of  bits  of  the  phase  excitations  is  not 
important  because  the  optimal  solution  implies  zero  phases  impressed  on  the  dipole  antennas). 
This  solution  for  digitized  excitations  is  derived  from  the  solution  for  "continuous"  excitations  by 
taking  the  digitized  excitations  that  are  closest  to  the  "continuous  excitations".  Fig.  11  shows  the 
3D-plot  of  the  synthesized  radiation  pattern  with  digitized  excitations  and  Figs.  12-14  shows  the 
total  radiation  field  in  the  H-,  D-  and  £-pIanes  for  both  the  digitized  and  the  "continuous" 
excitations. 


GENERATED  |E|  FOR 
DIGITIZED  EXCITATION 

SM(0)cos(«)  siN<e)siN(;0) 


Prescribed  beam  width  =  30° 
Number  of  elements  =  81 


The  synthesized  total  radiation  pattern  corresponding  to  grid  I  for  the  case  of  2-bits 
amplitude  excitations. 


RADIATION  PATTERN 
IN  H-PLANE 


=  PRESCRIBED  |  E  | 

=  SYNTHESIZED  |  E  | 

•  =  GENERATED  |  E  |  FOR 
DIGITIZED  EXCITATION 


SIN(0) 

Prescribed  beam  width  =  30“ 
Number  of  elements  =  81 


The  total  radiation  pattern  in  //-plane  of  the  prescribed,  digitized  and  "continuous" 
synthesized  radiation  field  coirespionding  to  grid  I  for  tlie  case  of  2-bits  amplitude 
excitations. 
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RADIATION  PATTERN 
IN  D-PLANE 


Fig.  13; 


—  =  PRESCRIBED  |  E  | 

=  SYNTHESIZED  |  E  | 

-  GENERATED | E | FOR 

DIGITIZED  EXCITATION 


Prescribed  beam  width  =  30° 
Number  of  elements  =  81 


The  total  radiation  pattern  in  Z)-p!ane  of  the  prescribed,  digitized  and  "continuous" 
synthesized  radiation  field  corresponding  to  grid  I  for  the  case  of  2-bits  amplitude 
excitations. 


t 

LU 


RADIATION  PATTERN 
IN  E-PLANE 


—  =  PRESCRIBED  |  E  | 

=  SYNTHESIZED  |  E  | 

--  =  GENERATED  |  E  |  FOR 
DIGITIZED  EXCITATION 


Number  of  elements  =  81 


Fig.  14:  The  total  radiation  pattern  in  F-planc  of  the  prescribed,  digitized  and  "continuous" 

synthesized  radiation  field  corresponding  to  grid  I  for  the  case  of  2-bits  amplitude 
excitations. 
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From  these  plots  we  see  that  for  the  case  of  the  digitized  amplitude  excitations  the  approximation 
of  the  main  beam  is  nearly  as  accurate  as  for  the  ccatinuous  case  while  the  sidelobes  of  the 
radiation  pattern  for  the  digitized  amplitude  excitations  are  higher  than  the  ones  for  the 
continuous  case.  For  the  digitized  case  we  find  NERR=42%  which  is  slightly  higher  than  the 
value  of  NERR  for  the  continuous  case. 
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6  CONCLUSION 

A  novel  method  of  array  antenna  pattern  synthesis  is  presented  that  is  based  on  an  optimization 
technique.  The  complex  excitation  values  of  the  current  distributions  of  the  array  elements  are 
taken  as  parameters.  The  radiated  far-field  pattern  is  matched  to  one  that  is  prescribed  by 
minimizing  a  global  mean-square  error.  The  method  is  applicable  to  arrays  with  unequaUy- 
spaced,  non-identical  elements  (for  which  a  factorization  of  the  radiation  pattern  into  an  element 
factor  and  an  array  factor  is  impossible).  Numerical  experiments  show  the  feasibility  of  the 
method. 
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DERIVATION  OF  CONDITION  FOR  OPTIMUM  EXCITATION 

In  this  appendix  the  condition  for  the  optimum  excitations  is  derived  from 

Equation  (27).  To  this  end,  the  variational  analysis  presented  in  [31]  is  applied.  We  put 

c((0  j)  =  c°PK£o;/)  +  Bc((oJ)  for  1=1,. (Al) 
where  5c(co;/)  is  an  arbitrary  deviation  from  Substitution  of  (Al)  into  (27)  leads  to 

N 

]|£“-£0»||„2=  ||£D»>||^2-2Re{  I  Jc°P‘*((o;f)+5c*(CO;/)]<£0~(|,to),e“(|,to;/)>) 

N  N 

+  I  ^  i  ^[cop'(co;/)+5c(a);/)][c«P'*((i);y)+&:*((i);/)J  <e~(§,co;/),e“(^,to;y)> 

’  N 

=  ERR2(£“0P')  +  2Re  I  dc\(aj) 

L  /=1 

N 

X  c°P'(a);/)<e“(|,(o;/),e~(|,co;y)>  - 
N  N 

+  I_^5c(0)/)5c*(co;y)<e“’(|.0);/),e~(^,0);y)>.  (A2) 

For  c°P'(co;/),  1=1,.. .,N,  to  be  a  minimum  of  ERR  (i.e.,  we  have  a  minimum  for  8c{(Oj)=0, 
/=1,...,A0,  the  terms  in  (A2)  that  are  linear  in  5c*(03  j)  have  to  compensate  each  other.  This  leads 
to  the  condition 

;^^c«P*(co;/)<e'”(|,co;/),e“(|,aj;y)>  =  <£O-(4,t0),e“(^,0);y)> 

for/=l,...7V.  (A3) 


TNO  report 


Appendix  B 


Page 

B.l 


DERIVATION  OF  EXPLICIT  EXPRESSION  FOR  Fp^ 

In  this  appendix  the  explicit  expression  for  Fp^(x^,(£))  in  the  case  of  a  constant  weighting  function 
w(4)=l  is  derived.  From  (31)  and  (32)  we  then  have 


Fp^ix^M  = 


(Bl) 


in  which  rj=A;j->j.  Integration  of  the  first  term  T1  in  (Bl)  gives 


Tl  =  (a)Po)25, 


exp(Jko^s^s)<^  -  47t(0)|ao)28  -sin(I:oR)/(*oR), 


(B2) 


where  /?=lrl.  The  integration  of  the  second  term  T2  is  found  by  using 


T2  =  -(copo)‘ 


§€Q 


4p4,expO'I:o^/.)dQ  =  [(topo)^/^'o^]3p3. 


exp(jk(^/,)(ia,  (B3) 


where  9_  denotes  partial  differentiation  with  respect  to  r  Using  the  results  of  (B2)  into  (B3)  leads 


to 


T2  =  4Kl(0)PfiMo^dpd^[sm(koRmoR)] 

-  47t[(copo)^/(^-o^)^]  |5p,[cos(^(A)-sin(I:o/?)/(^R)] 

-  {3[cos(i(/?)  -  sin(I:o/?)/(^oR)]  +  (I:o/?)sin(I:of?)}  1.  (B4) 


Summing  the  two  terms  T1  and  T2  gives  for  (Bl) 
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Fp^ixjM  =  4ji(co^o)2  I  \sin(koR)KkoR) 


+  (koR)-^[cos(koR)-smikMkoR)] 


For  the  case  r=0  {x=y)  the  expression  for  can  be  found  by  taking  the  limit  R—^  in  (B5).  A 
more  straightforeward  manner  is  the  substitution  of  r=0  in  formula  (Bl)  which  entails  that  the 
exponential  function  becomes  one,  and  hence 


Fp,(X,X.(0)  =  (COPo)2(5p,-^^)da 


For  the  determination  of  (B6)  we  introduce  spherical  coordinates  0  and  0,  where  0  is  the  angle 
between  ^  and  the  positive  Xj-axis  and  0  the  angle  that  the  projection  of  ^  on  the 
makes  with  the  positive  X;-axis.  In  terms  of  these  coordinates,  the  unit  vector  in  the  direction  of 
observation  ^  is 

^  =  [sin(0)cos(0),sin(0)sin(0),cos(0)]’*’,  (B7) 

and  the  Jacobian  is  sin(0).  Substitution  of  (B7)  into  (B6)  gives  for  the  expressions  of  (B6)  for 
which  p=q\ 


F/;(xjr,co)  =  (copo)2j^_^  d0  J  [l-sin2(0)cos2(0)]sin(0)d0=  (8/3)7t(a)Mo)2,  (B8) 


F22(x jr.co)  =  (0)p<,)2 J  d0 


and 


n 

[l-sin2(0)sin2(0)jsin(0)d0=  (8/3)7r(  topo)2, 

0=0 
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r 

2n  K 

Fjj(xjc,(o)  =  do  [l-cos2(0)]sin(0)d6  =  (8/3)7C(coHoF 

J  0=0  J  0=0 


(BIO) 


For  the  case  that  p  ^  ^  the  components  of  the  expression  in  (B6)  are  zero.  Combining  the  results 
(B8)-(B10)  we  find 


F.,(xj-,(0)  =  (8/3)Jt((4io)28... 


(BID 
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ANALYTICAL  EVALUATION  OF  INTEGRAL 

In  this  Appendix  we  perform  the  analytical  evaluation  of  integral  (42) 


r 


INTI>(y),co]  =  2(copo)2 


|6Q' 


^3  {  //[4j-cos(5)]  )  exp[-7l:o^^^(/)]cin, 


for7=l,...A'£>. 


(Cl) 


For  the  determination  of  (Cl)  we  introduce  spherical  coordinates  0  and  0,  where  6  is  the  angle 
between  ^  and  the  positive  Xj-axis  and  0  the  angle  that  the  projection  of  ^  on  the  XjrX2-plane 
makes  with  the  vector  y{J).  This  choice  with  respect  to  the  angle  0  leads  to  a  simple  expression 
for  the  inner  product  in  the  exponent  of  the  exponential  function.  In  terms  of  these  coordinates, 
the  unit  vector  in  the  direction  of  observation  ^  is 


^  =  [sin(0)cos(o),sin(0)sin(0),cos(0)p.  (C2) 

the  Jacobian  is  sin(0)  and  for  the  inner  product  we  find 

^,_Vj(/)=ly(/)lsin(0)cos(0).  (C3) 

Substitution  of  (C2)  and  (C3)  into  (Cl)  gives 


INTL>-(7).co]  =  2(to^o)2 


2n 

0=0 


d0 


ic/2 

0=0 


cos(0)  j//[cos(0)-cos(5)] ) 


exp[-yj:oly(/)lsin(0)cos(0)]sin(0)d0,  for  7=1  ....JVD. 


(C4) 


Equation  (C4)  is  rewritten  as 


INT[>-(7).co]  =  4((opo)2 


rs 

0=0 


sin(0)cos(0)d0 


11 

0=0 


exp[-y^oly(7)lsin(0)cos(0)]d0. 


(C5) 
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Using  the  theory  of  Bessel  functions  (see  [32])  the  integration  with  respect  to  0  leads  to 


INT[y(/),(o]  =  47i((Diio)2j  ^  ^sin(6)cos(9>/o[<:Qly(J)lsin(6)]d6, 
for  7=1,. ..//D. 


in  which  Jq  is  the  Bessel  function  of  the  first  kind  and  of  order  zero.  Next,  using  the  property  (see 
[32]) 


Uo{t)  =  -  [t/;(0] 

At 


we  find 


47l(copo)2Q/;(e) 

INT|>(7),(o]  =  - —  .  for  7=1,...7VD 

(^■0l3'(-OI)2 


where  Jj  is  the  Bessel  function  of  the  first  kind  of  order  one  and  Q=it(jly(/)lsin(5). 
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